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In brief 

Human antibody responses to SARS- 
CoV-2 differ between vaccination and 
infection, with mRNA vaccination 
inducing more productive lymph node GC 
responses and several vaccine types 
stimulating IgG antibodies capable of 
recognizing a broader range of viral 
variants. 
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SUMMARY 


During the SARS-CoV-2 pandemic, novel and traditional vaccine strategies have been deployed globally. We 
investigated whether antibodies stimulated by mRNA vaccination (BNT162b2), including third-dose boosting, 
differ from those generated by infection or adenoviral (ChAdOx1-S and Gam-COVID-Vac) or inactivated viral 
(BBIBP-CorV) vaccines. We analyzed human lymph nodes after infection or mRNA vaccination for correlates 
of serological differences. Antibody breadth against viral variants is lower after infection compared with all 
vaccines evaluated but improves over several months. Viral variant infection elicits variant-specific anti- 
bodies, but prior MRNA vaccination imprints serological responses toward Wuhan-Hu-1 rather than variant 
antigens. In contrast to disrupted germinal centers (GCs) in lymph nodes during infection, MRNA vaccination 
stimulates robust GCs containing vaccine mRNA and spike antigen up to 8 weeks postvaccination in some 
cases. SARS-CoV-2 antibody specificity, breadth, and maturation are affected by imprinting from exposure 
history and distinct histological and antigenic contexts in infection compared with vaccination. 





INTRODUCTION 


The urgent need for countermeasures against the coronavirus 
disease 2019 (COVID-19) pandemic has spurred the rapid devel- 
opment of severe acute respiratory syndrome coronavirus 2 
(SARS-CoV-2) vaccines of diverse formulations. mRNA vaccines 
BNT162b2 (BioNTech-Pfizer) and mRNA-1273 (Moderna/NIAID) 


have demonstrated high efficacy and safety in clinical trials for 
COVID-19 prevention (Baden et al., 2021; Polack et al., 2020; 
Walsh et al., 2020). Additional COVID-19 vaccines including 
adenoviral vectored vaccines ChAdOx1-S (Astra Zeneca) (Voy- 
sey et al., 2021), Ad26.COV2.S (Johnson & Johnson) (Sadoff 
et al., 2021), and Gam-COVID-Vac (Sputnik V) and inactivated 
viral vaccines such as BBIBP-CorV (Sinopharm) also have 
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reported efficacy. Correlates of vaccine-elicited protection from 
COVID-19 are the titers of neutralizing antibodies to SARS- 
CoV-2, and the concentration of antibodies binding to spike or 
receptor-binding domain (RBD) (Earle et al., 2021; Gilbert 
et al., 2022; Khoury et al., 2021; Róltgen and Boyd, 2021). 
Most neutralizing antibodies target the RBD and prevent binding 
to angiotensin-converting enzyme 2 (ACE2) receptor (Greaney 
et al., 2021a; Yuan et al., 2021). Current SARS-CoV-2 vaccines 
all contain or induce the expression of antigens similar to those 
of the early Wuhan-Hu-1 viral isolate, but differ in elicited binding 
and neutralizing antibody responses, with higher responses from 
mRNA vaccines compared with adenovirus-vectored or inacti- 
vated virus vaccines (Dashdorj et al., 2021a, 2021b). It remains 
to be determined precisely how the immune system responds 
to mRNA and other vaccine platforms compared with SARS- 
CoV-2 infection. Data from RBD variant antigen yeast display 
and pseudotyped virus neutralization show that RBD epitope 
targeting by polyclonal serum antibodies is narrower in infected 
patients compared with mRNA-1273 vaccinees (Greaney 
et al., 20210). 

Several SARS-CoV-2 variants of concern with mutations in the 
spike gene have emerged and spread globally, with differing 
abilities to evade neutralizing antibody responses elicited by Wu- 
han-Hu-1 infection or vaccination. The most immune-evasive 
variants, including the recent Omicron variant, have alterations 
in epitopes containing amino acid E484 (Garcia-Beltran et al., 
2021; Greaney et al., 2021a; Hoffmann et al., 2021). The appear- 
ance of virus variants, waning antibody levels after infection or 
vaccination (Falsey et al., 2021; Levin et al., 2021), and break- 
through infections in previously immunized individuals (Keehner 
et al., 2021) indicate that periodic vaccine boosting of immunity 
to SARS-CoV-2 is warranted. Third doses of mRNA-1273 (Chu 
et al., 2021) and BNT162b2 (Falsey et al., 2021) administered 
several months after the second dose prompt an increase in 
neutralizing antibodies greater than the peak following initial 
vaccination doses. mRNA-1273 vaccination followed by mRNA 
booster vaccines expressing Beta spike gives higher neutralizing 
titers to Wuhan-Hu-1-like SARS-CoV-2 compared with the Beta 
variant (Choi et al., 2021; Wu et al., 2021), suggesting that some 
degree of immune imprinting, or preferential responses to the 
viral variants initially encountered by the immune system, may 
affect the development of antibodies against new viral variants 
(Wheatley et al., 2021). 

Germinal center (GC) responses in human lymphoid tissues 
enable antibody affinity maturation and durable serological and 
memory B cell responses although extrafollicular B cell re- 
sponses are also reported (Elsner and Shlomchik, 2020; Lam 
et al., 2020; Woodruff et al., 2020). The degree to which SARS- 
CoV-2 infections or different vaccines stimulate GC responses 
and differ in factors such as the quantity, persistence, and local- 
ization of antigen in lymph nodes (LNs) and other lymphoid tis- 
sues are important open questions. Approaches such as fine- 
needle aspiration (FNA) are being increasingly used to study 
LN-derived cells from healthy human subjects (Havenar-Daugh- 
ton et al., 2020; Lederer et al., 2022; Turner et al., 2021). Disrup- 
ted LN GCs have been reported in autopsies of deceased 
COVID-19 patients (Haslbauer et al., 2021; Kaneko et al., 
2020), while elevated frequencies of GC B cells are seen after 
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mRNA vaccination in healthy individuals (Turner et al., 2021) 
and lower GC B cell frequencies after mRNA vaccination of 
immunocompromised individuals (Lederer et al., 2022). To 
date, no direct comparison of LN GC histology and cellular 
composition, combined with measurement of viral or vaccine an- 
tigen quantity, persistence, and distribution in draining LN sites 
of COVID-19 patients and vaccinees has been reported. 

Here, we compare antibody responses in BNT162b2 mRNA 
vaccine recipients following first, second, and third vaccine 
doses to antibody responses of COVID-19 patients. We find dif- 
ferences in the magnitude, isotype profiles, SARS-CoV-2 spike 
domain specificity, and breadth of binding antibody responses 
to a panel of nine viral variants in addition to Wuhan-Hu-1. 
Anti-RBD IgG binding to SARS-CoV-2 variants of concern and 
interest in recipients of four different vaccines (BNT162b2, ChA- 
dOx1-S, Gam-COVID-Vac, and BBIBP-CorV) and in COVID-19 
patients shows greater binding breadth for viral variant RBDs 
following all vaccines compared with Wuhan-Hu-1 infection. 
We quantify a strong imprinting effect of prior vaccination with 
Wuhan-Hu-1 spike antigen on antibody specificities following 
breakthrough infection with viral variants. Histological analysis 
of draining LN shows marked impairment of GCs in severe 
COVID-19 compared with mRNA vaccination, higher quantities, 
and persistence of spike antigen accumulated in the GCs of 
mRNA vaccinees and detectable vaccine RNA in GCs for up to 
2 months post-second dose. 


RESULTS 


Magnitude and waning of anti-SARS-CoV-2 IgG 

following BNT162b2 vaccination and response to third- 
dose boost 

We measured anti-SARS-CoV-2 antibodies for nucleocapsid 
(N), full spike, and RBD in Stanford BNT162b2 study participant 
plasma samples using multiplexed electrochemiluminescence 
(ECL) assays (Meso Scale Discovery, MSD), in WHO binding 
antibody units (BAU). The first and second vaccine doses 
were at days 0 and 21, with the third-dose boosting at approx- 
imately 9 months. Plasma samples were collected in a time 
course up to 7 months after the first dose and up to 1 month 
after the third dose. Four of the 59 vaccine recipients had a his- 
tory of SARS-CoV-2 reverse-transcription quantitative polymer- 
ase chain reaction (RT-qPCR)-confirmed infection (CoV-2+) 
prior to vaccination. IgG for spike protein and RBD in vaccinees 
reached their initial peak at day 28 after the first dose (Figures 
1A and 1B). IgG binding to spike and RBD was highly corre- 
lated with SARS-CoV-2 neutralization titers (Arunachalam 
et al., 2021) (Figure 1C). By 9 months, spike-specific IgG had 
decreased approximately 20-fold from the maximum, but the 
third-dose boost raised IgG concentrations above the prior 
peak within 1 week. IgG specific for N protein, which is not en- 
coded in the vaccine, was negative throughout the study in 54 
of the previously uninfected vaccinees, but one participant se- 
roconverted for anti-N IgG between days 90 and 210 after the 
prime, indicating a breakthrough infection (Figure 1A). CoV-2+ 
vaccinees had accelerated RBD and spike IgG responses after 
the first dose, and detectable anti-N IgG was unaffected by 
vaccination (Figure 51А). 
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Figure 1. Magnitude and duration of anti-SARS-CoV-2 IgG following BNT162b2 vaccination and third-dose boost 

(A) Anti-SARS-CoV-2 N, RBD, and spike (S) antibody responses are shown for plasma samples from individuals who received BNT162b2 prime (DO, n = 59 
individuals), second dose (D21, n = 58 individuals), and third dose (around month 9, n = 36 individuals) vaccination. Box-whisker plots of the anti-SARS-CoV-2 IgG 
concentrations in WHO BAU/mL show the interquartile range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range 
below the 2596 quantile and above the 7596 quantile. Red dashed lines indicate the cutoff values for the positivity of each assay (MSD, package insert). 

(B) Heatmap showing the development of antibody responses in longitudinal samples collected at DO, D7, D21, D28, D42, and D90/120 time points post-prime 
vaccination (x axis). WHO BAU/mL Ig concentrations are displayed for study participants sorted by age (y axis, color coded). Rows are labeled on the right with 
“CoV-2+” for participants with a previous SARS-CoV-2 RT-qPCR positive test result. 

(C) Correlations between anti-RBD and anti-spike IgG binding antibody concentrations in WHO BAU/mL and SARS-CoV-2 virus neutralization assays. Spearman 
rank correlation (coefficient = Rho, displayed in the plot for each assay comparison) was used to assess the strength of correlation between binding antibody 
concentrations and virus neutralization results. Red dashed lines indicate the cutoff values for the positivity of each assay (MSD, package insert). 
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BNT162b2 recipients had weak IgM and IgA responses to 
spike and RBD compared with their IgG responses. Robust IgG 
responses were seen in all age groups (Figures 1B and S1A- 
S1C). Convalescent COVID-19 patients and BNT162b2 vacci- 
nees had similar low saliva IgG concentrations for spike and 
RBD, several orders of magnitude lower than those detected in 
plasma (Figure S1D). As in plasma, saliva IgG peaked at 1 week 
after third-dose boosting, at higher values than the peak after 
the second dose (Figure S1D). Reported side effects after vacci- 
nation showed no relationship to plasma IgG responses (Figures 
S2A and S2B). 


BNT162b2 vaccination and Wuhan-Hu-1 SARS-CoV-2 
infection stimulate distinct antibody isotypes and 
endemic coronavirus antibody responses 

Severe COVID-19 stimulates higher SARS-CoV-2-specific anti- 
body titers than asymptomatic infection or mild illness (Long 
et al., 2020; Roltgen et al., 2020). We compared antibody isotype 
concentrations specific for spike and RBD in COVID-19 patients 
(Stanford cohort 1 of this study) from the initial months of the 
pandemic (Róltgen et al., 2020) with the responses of the Stanford 
BNT162b2 vaccinees (Figures 2A and 2B). Patients were classified 
as outpatients, admitted patients not requiring care in the intensive 
care unit (ICU), ICU patients, and those who died from their illness. 
Stanford BNT162b2 vaccinee RBD and spike IgG concentrations 
were comparable to those of severely ill patients and higher than 
those of mildly or moderately ill patients for anti-RBD antibodies 
at day 42 (Figures 2A and 2B). The BNT162b2 vaccine induced a 
highly IgG-polarized serological response with minimal IgM- and 
IgA-binding spike and RBD (Figures 2A and 2B). Principal compo- 
nent analysis (PCA) showed clustered and homogeneous SARS- 
CoV-2 spike and spike-domain-specific serological responses in 
BNT162b2 vaccinees compared with infected patients, as quanti- 
fied by smaller distances for vaccinated participants from the 
group centroid (Figures 2C and 2D). 

SARS-CoV-2 vaccinees and COVID-19 patients showed 
boosting of SARS-CoV-1 spike antibodies, but infected patients 
showed greater boosting of spike IgG and IgA for endemic hu- 
man betacoronaviruses OC43 and HKU1 (Figures S3A and 
S3B). The third BNT162b2 vaccine dose further increased 
vaccinee titers to SARS-CoV-1, OC43, апа HKU1 (Figure 530). 
Antibodies to the spike antigens of the endemic human alpha- 
coronaviruses NL63 and 229E were not boosted (Figure S3). 


с? CelPress 


N ACC 


Greater breadth of IgG binding to viral variants following 
BNT162b2 vaccination compared with infection with 
Wuhan-Hu-1 SARS-CoV-2 

Immune-evasive SARS-CoV-2 variants have spread globally 
(Harvey et al., 2021; Plante et al., 2021; Róltgen and Boyd, 
2021). We compared plasma IgG responses to the RBDs of 
nine different SARS-CoV-2 variants of concern and interest in 
BNT162b2 vaccinees and COVID-19 patients, using multiplexed 
MSD ECL assays. For RBD antigens from Epsilon, Kappa, 
B.1.526.2, B.1.214.2, Alpha, Eta/lota, Gamma, P.3, and Beta 
variants, both vaccinee and infected patient IgG showed the 
greatest decrease in binding to Beta, Gamma, and P.3 variants 
relative to Wuhan-Hu-1 (Figure 3A). To quantify the differences 
in variant RBD binding by vaccinee and patient plasma IgG, 
we calculated the ratios of anti-RBD IgG concentrations for Wu- 
han-Hu-1 compared with viral variants, with higher ratios indi- 
cating greater binding of Wuhan-Hu-1 RBD compared with 
variant RBD (Figure 3B). COVID-19 patients showed a greater 
IgG binding bias for Wuhan-Hu-1 RBD compared with variant 
RBDs in the initial weeks post-onset of symptoms; in contrast, 
BNT162b2 vaccinee IgG had relatively greater breadth of bind- 
ing to variant RBDs and less preference for Wuhan-Hu-1 RBD. 
Over time, infected patient plasma samples showed improve- 
ment in variant RBD binding relative to Wuhan-Hu-1 RBD, sug- 
gesting evolution of the antibody response through at least 
7 weeks post-onset of symptoms (Figure 3). BNT162b2 
vaccinee IgG Wuhan-Hu-1 to variant RBD binding ratios did 
not change from day 21 onward. The greater breadth of variant 
RBD binding (including the Delta variant) by vaccinee IgG 
compared with COVID-19 patient IgG was seen in a second, in- 
dependent validation cohort (Stanford cohort 2) of predomi- 
nantly mildly ill COVID-19 patients. Greater Wuhan-Hu-1 to 
variant RBD IgG binding ratios were found in weeks 2-3, month 
1, month 3, and month 7 in infected patients compared with vac- 
cinees (Figures S4A and S4B), with improvement in variant 
recognition over time in the infected patients. Notably, the 
increased breadth of vaccinee IgG compared with COVID-19 
patient IgG binding to viral variant antigens was greatest for 
RBD, the main target of neutralizing antibodies, and was 
decreased or not detected when whole spike antigens were 
tested (Figure S4C). Functional blocking of ACE2 binding to 
RBD was concordant with the RBD-specific IgG concentrations 
measured in these populations (Figures S4A and S4D). 


Figure 2. BNT162b2 vaccination and SARS-CoV-2 infection elicit distinct antibody profiles 





(A and B) Anti-SARS-CoV-2 N, RBD, and spike (S) IgM, IgG, and IgA antibody responses are shown for individuals who received BNT162b2 prime (DO) and 
second (D21) vaccination doses and for COVID-19 patients. (A) The heatmap shows the development of antibody responses in longitudinal samples from 
vaccinees/patients collected at DO, D7/week 1, D21/weeks 2 and 3, D28/week 4, D42/weeks 5 and 6, and D90/120/- week 7 after vaccination/COVID-19 
symptom onset (x axis). The color scale encodes the median values of Іод:о WHO BAU/mL Ig concentrations. (B) Box-whisker plots show the development of 
antibody responses in longitudinal samples from vaccinees/patients collected at DO, D7/week 1, D21/weeks 2 and 3, D28/week 4, D42/weeks 5 and 6, and D90/ 
120/> week 7 after vaccination/COVID-19 symptom onset (x axis). Box-whisker plots show the interquartile range as the box and the whisker ends as the most 
extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 quantile. Statistical test: pairwise Wilcoxon rank sum test with 
Bonferroni correction. *p < 0.05, **р < 0.01, ***p < 0.001. Individuals were classified as outpatients (Outpt) and hospital-admitted patients (Admit); ICU patients 
and those who died from their illness (Death); and vaccinees who had (CoV-2+) or had not had a positive SARS-CoV-2 test in the past. 

(C) PCA of anti-SARS-CoV-2 RBD, N-terminal domain, and S (but not N) IgM, IgG, and IgA concentrations across BNT162b2 vaccinees and Wuhan-Hu-1-in- 
fected Stanford COVID-19 patient cohort 1 at different time points after vaccination/COVID-19 symptom onset visualized on a consistent PCA reference created 
using D21/weeks 2 and 3 as a reference time point. 

(D) Distribution of Euclidean distances between BNT162b2 vaccinee samples and their centroid, compared with Wuhan-Hu-1-infected Stanford COVID-19 
patient cohort 1 samples and their centroid, at different time points after vaccination/COVID-19 symptom onset. 
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Figure 3. Greater breadth of IgG binding to SARS-CoV-2 RBD variants following BNT162b2 vaccination compared with infection with Wuhan- 
Hu-1 SARS-CoV-2 

Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for Stanford individuals who received BNT162b2 vaccination and for Wuhan-Hu- 
1-infected COVID-19 Stanford patient cohort 1 at different time points after vaccination/COVID-19 symptom onset. Box-whisker plots show the interquartile 
range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 quantile. 
Significance between patient and vaccinee groups were tested with two-sided Wilcoxon rank sum test. *p < 0.05, **р < 0.01, ***p < 0.001 

(A) Anti-RBD IgG concentrations. 

(B) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration. 


Improved IgG binding to viral variants is consistent across 
four COVID-19 vaccines (BTN162b2, ChAdOx1-S, Gam- 
COVID-Vac, and BBIBP-CorV) compared with infection 
Several COVID-19 vaccines, including mRNA, viral vector- 
based, and inactivated virus vaccines, have been approved 
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for use internationally. Varying efficacy and antibody re- 
sponses from the vaccines have been reported (Baden 
et al., 2021; Dashdorj et al., 2021a, 2021b; Polack et al., 
2020; Sadoff et al., 2021; Voysey et al., 2021). We compared 
IgG responses in Stanford COVID-19 cohort 2 patients and 
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Figure 4. Greater breadth of IgG binding to SARS-CoV-2 variant RBDs following vaccination with four different vaccines compared with 
infection with Wuhan-Hu-1 SARS-CoV-2 

Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for individuals who received BNT162b2 (BioNTech-Pfizer), ChAdOx1-S (Astra 
Zeneca), Gam-COVID-Vac (Sputnik V), and BBIBP-CorV (Sinopharm) vaccination and for Wuhan-Hu-1-infected COVID-19 Stanford patient cohort 2 within 
1 month and around 3 months after vaccination/COVID-19 symptom onset. Box-whisker plots show the interquartile range as the box and the whisker ends as the 
most extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 quantile. Significance between groups was tested with 


pairwise Wilcoxon rank sum test with Bonferroni correction. *р < 0.05, **p < 0.01, ***р < 0.001 


(A) Anti-RBD IgG concentrations. 
(B) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration. 


BNT162b2 vaccinees with those of participants in a Mongo- 
lian observational study deploying four different COVID-19 
vaccines: the MRNA vaccine BNT162b2, adenoviral vectored 
vaccines ChAdOx1-S (Astra Zeneca) and Gam-COVID-Vac 
(Sputnik V), and an alum-adjuvanted, inactivated viral vaccine 
BBIBP-CorV (Sinopharm). RBD-specific IgG concentrations 
for Wuhan-Hu-1 and all viral variants measured (Epsilon, 
Kappa, B.1.526.2, Delta, Alpha, Eta/lota, Gamma, Р.З, 
and Beta) differed greatly between vaccine groups, with 
BNT162b2 eliciting the highest antibody levels, followed by 
Astra Zeneca, Sputnik V, and Sinopharm vaccination (Fig- 
ure 4A). IgG concentration differences between vaccines 
were significant for most viral variant RBDs. Stanford 
BNT162b2 vaccinees compared with Mongolian BNT162b2 


vaccinees had higher IgG concentrations at early time points, 
likely due to differences in the timing of sample collection 
(Stanford days 28 and 90; Mongolian participants variable 
time points before 1 and 3 months) (Figure 4A). Despite the 
different vaccine compositions and magnitudes of antibody 
responses, all four vaccines elicited IgG with relatively greater 
breadth of viral variant RBD binding compared with that of in- 
fected patients (Figure 4B). 


Variant-specific serological responses following Alpha 
and Delta SARS-CoV-2 infection and immune imprinting 
after vaccination 

Immune imprinting, a phenomenon in which primary exposure to 
an antigen forms epitope-specific B cell memory and affects 
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Figure 5. Variant-specific serological signature following Alpha and 
Delta SARS-CoV-2 infection 

(A) Anti-Wuhan-Hu-1 to variant RBD IgG concentration ratios are shown for 
individuals with primary SARS-CoV-2 Alpha or Delta variant infection (upper 
panels) or secondary variant infection after vaccination (lower panels). Box- 
whisker plots show the interquartile range as the box and the whisker ends as 
the most extreme values within 1.5 times the interquartile range below the 25% 
quantile and above the 75% quantile. 

(B) Anti-SARS-CoV-2 variant IgG binding preference levels of BNT162b2 
vaccinees on day 28 postvaccination and of previously vaccinated or non- 
vaccinated individuals infected with the SARS-CoV-2 Delta variant. 


future B cell and antibody responses against variant epitopes, 
has been studied in influenza infection and vaccination. 
COVID-19 patients and BNT162b2 vaccinees who were only 
exposed to Wuhan-Hu-1 antigens in this study exhibit a consis- 
tent hierarchy in IgG binding concentrations to the different 
SARS-CoV-2 variant RBDs relative to the Wuhan-Hu-1 RBD, 
decreasing from Epsilon, Kappa, B.1.526.2, Delta, Eta/lota, 
Gamma, P.3 to Beta (Figure S5A). To test for imprinting of the 
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serological response to variant RBDs, we first analyzed the 
ratios of Wuhan-Hu-1 to variant RBD IgG concentrations in 
COVID-19 patients who were infected with Alpha or Delta vari- 
ants, confirmed by allele-specific RT-qPCR testing or viral 
sequencing. IgG from Alpha- or Delta-variant-infected patients 
with no history of COVID-19 vaccination or prior SARS-CoV-2 
infection preferentially bound Alpha and Delta variant RBDs, 
respectively, compared with Wuhan-Hu-1 RBD (Figure 5A, up- 
per panels). Delta infection also elicited higher IgG concentra- 
tions to other variant RBDs containing L452R such as Epsilon 
and Kappa, compared with Wuhan-Hu-1 (Figure 5A, upper right). 
PCA of variant RBD-specific IgG responses in vaccinees and 
variant-infected patients (Figure S5B) highlights the distinct 
serological responses elicited by infection with the variant 
viruses. To test whether prior exposure to one SARS-CoV-2 
RBD variant causes imprinting of humoral immunity, we 
analyzed plasma from individuals vaccinated with Wuhan-Hu- 
1-like antigens and subsequently infected with Alpha or Delta 
variants (Figure 5A, lower panels). Despite breakthrough infec- 
tion with Alpha or Delta viral variants, the vaccinated individuals 
showed patterns of IgG binding to viral variant RBDs similar to 
those of individuals exposed to only Wuhan-Hu-1. We quantified 
the degree of imprinting of IgG specificity by log transforming the 
ratios of IgG binding to pairs of antigens (e.g., Wuhan-Hu-1 RBD 
compared with Delta RBD) for individual samples, then rescaling 
to range from —100% to +100% corresponding to the maximal 
preference for each antigen observed in other plasma speci- 
mens, including those from individuals exposed only to a single 
antigen variant (Figure 5B). 


LN GC impairment in severe SARS-CoV-2 infection but 
robust development following SARS-CoV-2 mRNA 
vaccination 

The differences in viral variant RBD IgG binding between SARS- 
CoV-2-infected patients and recipients of the four COVID-19 
vaccines suggest that the organization of the humoral immune 
responses in secondary lymphoid tissues may differ between 
infection and vaccination, potentially due to direct effects of 
the viral infection, differences in innate immune stimuli between 
vaccination and infection, or the quantity or localization of viral 
antigens, among other possibilities. Previous studies have re- 
vealed a loss of GCs and a reduction in BCL6* GC B cells in se- 
vere acute SARS-CoV-2 infection, raising the possibility that 
humoral responses may be altered or subverted by the virus (Ka- 
neko et al., 2020). It is unclear whether draining LN immune re- 
sponses to SARS-CoV-2 infection in the lungs differ from those 
elicited in axillary LNs following deltoid intramuscular mRNA 
vaccination. To compare GC architecture in response to 
SARS-CoV-2 infection and vaccination, we obtained peribron- 
chial LN tissues from seven COVID-19 patients and three control 
autopsy cases as well as axillary LN core needle biopsies of 
seven individuals vaccinated with mRNA-1273 or BNT162b2. 
Importantly, core needle biopsy sample tissue volumes were 
suitable for the assessment of LN histoarchitecture. Vaccinee 
axillary LN biopsies were from the ipsilateral (same-side) arm 
vaccinated. Controls were thoracic LNs from individuals who 
succumbed to prepandemic non-COVID-19 pneumonias and 
contain GCs likely due to ongoing adaptive immune responses 
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elicited by other antigens. LN histology for COVID-19 patients 
and vaccinees was evaluated with four-color codetection by in- 
dexing (CODEX) immunofluorescence analysis for CD20, CD3, 
BCL6, and CD21, which are markers of B cells, T cells, GC B 
cells (or T follicular helper [Tfh] cells), and follicular dendritic 
cells, respectively (Figures бА and 56А), as well as by single-co- 
lor immunohistochemical stains for these markers and the Tfh 
cell marker PD-1 (Figure 6B). GCs were poorly formed in the 
severely ill COVID-19 patient peribronchial LNs compared with 
the axillary LNs of vaccinees, with disrupted CD21* follicular 
dendritic cell networks and decreased BCL6* cells (including 
GC B cells and Tfh cells) and PD-1* cells (consistent with Tfh 
cells) within GCs (Figures 6A-6E). Disruption of CD21* follicular 
dendritic cell networks was seen in both primary and secondary 
follicles in LNs from COVID-19 patients (Figure S6B). mRNA 
vaccination was associated with follicular hyperplasia with fully 
developed GC architecture, including robust induction of GC B 
cells, Tfh cells, and extensive follicular dendritic cell networks 
(Figures 6A and 6B). 


Prolonged detection of vaccine mRNA in LN GCs and 
spike antigen in LN GCs and blood following SARS- 
CoV-2 mRNA vaccination 

The biodistribution, quantity, and persistence of vaccine mRNA 
and spike antigen after vaccination and viral antigens after 
SARS-CoV-2 infection are incompletely understood but are likely 
to be major determinants of immune responses. We performed 
in situ hybridization with control and SARS-CoV-2 vaccine 
mRNA-specific RNAScope probes in the core needle biopsies 
of the ipsilateral axillary LNs that were collected 7-60 days after 
the second dose of mRNA-1273 or BNT162b2 vaccination and 
detected vaccine mRNA collected in the GCs of LNs on days 
7, 16, and 37 postvaccination, with lower but still appreciable 
specific signal at day 60 (Figures 7A-7E). Only rare foci of vac- 
cine mRNA were seen outside of GCs. Axillary LN core needle bi- 
opsies of nonvaccinees (n = 3) and COVID-19 patient specimens 
were negative for vaccine probe hybridization. Immunohisto- 
chemical staining for spike antigen in mRNA-vaccinated patient 
LNs varied between individuals but showed abundant spike pro- 
tein in GCs 16 days post-second dose, with spike antigen still 
present as late as 60 days post-second dose. Spike antigen 
localized in a reticular pattern around the GC cells, similar to 
staining for follicular dendritic cell processes (Figure 7B). 
COVID-19 patient LNs showed lower quantities of spike antigen 
but a rare GC had positive staining (Figure 7F). Immunohisto- 
chemical staining for N antigen in peribronchial LN secondary 
and primary follicles of COVID-19 patients (Figures 7F-71) was 
positive in five of the seven patients, with a mean percentage 
of nucleocapsid-positive follicles of more than 2596. 

Spike protein was detected in the plasma of 9696 of the vacci- 
nees at days 1-2 (median spike concentration of 47 pg/mL) and 
in 6396 at day 7 (median spike concentration of 1.7 pg/mL) after 
the prime vaccine dose. In contrast, spike antigen detection after 
the vaccine boost on day 21 was reduced, with half of the study 
participants being positive on days 1-2 (median spike concen- 
tration of 1.2 pg/mL) and only one individual on day 7 post-boost 
(Figure 7J). We suspected that high concentrations of spike-spe- 
cific antibodies developed by vaccinees within the first 2— 
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3 weeks after the prime vaccine dose could impede detection 
of spike antigen by competing for spike-binding sites with the 
anti-spike reagent antibodies in the antigen assay. To test this 
hypothesis, we added different concentrations of recombinant 
spike protein to spike-negative vaccinee plasma samples 
collected on day 0 (n = 3) and day 28 (n = 3) after the prime vacci- 
nation. While the recombinant spike protein could readily be de- 
tected in day O plasma samples, only high concentrations of the 
antigen led to a positive signal when mixed with the day 28 sam- 
ples (Figure 7K). We then mixed a spike-positive plasma sample 
collected one day after vaccination with spike-negative plasma 
samples collected on days 0, 21, 22-23, and 28 (n = 4 each). 
Spike antigen detection levels were high in the mix of day 1 
and day 0 samples, decreased in the mix of day 1 and day 21, 
and day 1 and day 22-23 samples, and below the cutoff for pos- 
itive in the mix of day 1 and day 28 samples (Figure 7L). Together, 
our results are consistent with spike-specific antibodies blocking 
the detection of the antigen in antigen capture-based assays. 


DISCUSSION 


One of the positive developments amid the global calamity of the 
SARS-CoV-2 pandemic has been the rapid design, production, 
and deployment of a variety of vaccines, including remarkably 
effective mRNA vaccines encoding the viral spike (Baden 
et al., 2021; Polack et al., 2020). We find that BNT162b2 vaccina- 
tion produces IgG responses to spike and RBD at concentra- 
tions as high as those of severely ill COVID-19 patients and fol- 
lows a similar time course. Unlike infection, which stimulates 
robust but short-lived IgM and IgA responses, vaccination 
shows a pronounced bias for IgG production even at early time 
points. These responses were similar across the adult age range 
in our study. The relative absence of IgM and IgA responses sug- 
gests a potent effect of the vaccine formulation in driving early 
and extensive IgG class-switching, potentially as a result of the 
reported T helper type 1-polarized CD4* T cell responses stimu- 
lated by vaccine components (Lederer et al., 2020; Lindgren 
et al., 2017; Pardi et al., 2018). Our data demonstrate that 
vaccinee plasma and saliva spike and RBD-specific IgG concen- 
trations decrease from their peak values by approximately 20- 
fold by 9 months after primary vaccination but quickly exceed 
prior peak concentrations in 7-8 days after boosting with a third 
vaccine dose. 

Correlates of immunological protection from SARS-CoV-2 
infection following vaccination or prior infection are still 
under investigation. Analysis of Moderna mRNA-1273 and 
AstraZeneca ChAdOx1-S responses highlights the overall simi- 
larity of correlate of protection results for spike-binding antibody 
and neutralizing antibody assays (Feng et al., 2021; Gilbert et al., 
2022). We compared spike or RBD-binding antibody responses 
with Wuhan-Hu-1 SARS-CoV-2 neutralization data in BNT162b2 
vaccinees and confirmed the high correlation of these assay re- 
sults, supporting the interpretation that sensitive, precise, and 
validated commercial multiplexed antigen-binding assays with 
a wide dynamic range, such as the MSD ECL assays in this 
study, will be valuable in providing standardized correlates of 
protection data for vaccines as the pandemic continues. Partic- 
ularly in the context of viral variants, it will be important to 
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2 tients versus mRNA vaccinees 
Ф (А) Representative LN GC histology of COVID-19 
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о COVID-19, a control autopsy patient who died from 
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Figure 7. Localization of SARS-CoV-2 proteins and vaccine mRNA in LNs 
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(A) Representative LN GC after mRNA vaccination showing hematoxylin and eosin staining (upper left), four-color Codex staining (lower left), in situ hybridization 
of a SARS-CoV-2 mRNA vaccine-specific probe (upper right [lower magnification] and middle right [greater magnification], and immunohistochemical (IHC) 
staining for spike antigen (lower right). Vaccine mRNA probe hybridization was visualized by colorimetric development with Fast Red chromogen, and positive 
IHC staining for spike antigen was visualized as granular brown color from 3,3'-diaminobenzidine (DAB) reagent. 
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determine whether predictions of vulnerability to infection or se- 
vere disease can be improved by adding data from other immu- 
nological assays, including T cell measurements. 

Differences in B cell responses to SARS-CoV-2 infection and 
vaccination may be reflected in the binding breadth of antibodies 
to different SARS-CoV-2 variants. We find that plasma of individ- 
uals who received prime/boost BNT162b2 vaccination, as well 
as individuals who received adenoviral vectored (ChAdOx1-S 
or Gam-COVID-Vac) or inactivated virus (BBIBP-CorV) COVID- 
19 vaccines show consistent patterns of binding to variant 
RBDs with modest decreases compared with Wuhan-Hu-1 
RBD binding. In contrast, COVID-19 patients produce antibody 
responses with significantly greater Wuhan-Hu-1 RBD binding 
preference and lower breadth of variant RBD binding. These dif- 
ferences between vaccinee and COVID-19 patient IgG variant 
antigen binding were greatest for the RBD, the target of most 
neutralizing antibodies and were diminished when full spike an- 
tigen with its greater number of non-neutralizing epitopes was 
tested. These results, covering many clinically relevant viral 
variant antigens and several vaccine modalities, are consistent 
with findings for RBD-binding IgG in mRNA-1273 vaccinees 
compared with infected patients (Greaney et al., 2021b). 
Notably, COVID-19 patients with Alpha or Delta variant infec- 
tions display characteristic serological profiles specific to the 
RBD of the infecting variant, indicating that SARS-CoV-2 variant 
serotyping may be useful for epidemiological studies of popula- 
tions to determine exposure to circulating SARS-CoV-2 variants. 
Both vaccinees and COVID-19 patients exposed to Wuhan-Hu-1 
antigens show the greatest decreases in antibody binding to 
RBD variants harboring E484 alterations, including Beta and 
Gamma. Although susceptibility to infection by viral variants is 
common to both vaccinated and convalescent populations, 
particularly as antibody titers decrease over time (Israel et al., 
2022; Levin et al., 2021), our findings lead to the prediction that 
antibodies derived from infection may provide somewhat 
decreased protection against virus variants compared with com- 
parable concentrations of antibodies stimulated by vaccination. 

As additional variants of SARS-CoV-2 appear over time, indi- 
viduals will acquire distinct immunological histories depending 
on which vaccines they received and which viral variants in- 


Cell 


fected them. The idea that “imprinting” by a prior antigen expo- 
sure can shape, either positively or negatively, the response to a 
subsequent variant is well established in studies of influenza vi- 
ruses and has been implicated in birth-year differences in sus- 
ceptibility to particular avian influenza viruses (Gostic et al., 
2016). We find that prior vaccination with Wuhan-Hu-1-like anti- 
gens followed by infection with Alpha or Delta variants gives rise 
to plasma antibody responses with apparent Wuhan-Hu-1-spe- 
cific imprinting manifesting as relatively decreased responses to 
the variant virus epitopes, compared with unvaccinated patients 
infected with those variant viruses. While current booster vacci- 
nations are still based on the Wuhan-Hu-1-like antigens, vaccine 
manufacturers are in the process of evaluating updated vaccine- 
encoding sequences from one or more circulating variants. Initial 
results from the third-dose boosting with Beta-spike-encoding 
mRNA vaccines after prior second-dose mRNA-1273 vaccina- 
tion are consistent with our findings of significant imprinting of 
serological responses by the first antigen encountered (Choi 
et al., 2021; Chu et al., 2021), indicating that vaccine-derived 
imprinting affects subsequent antibody responses stimulated 
by vaccination as well as infection. The extent to which vaccine 
boosting or infection with different variants will effectively elicit 
antibody responses to new epitopes or rather increase re- 
sponses to the epitopes of antigens encountered previously, 
asinthe "original antigenic sin" phenomenon described for influ- 
enza virus infection and vaccination (Arevalo et al., 2020; Zhang 
et al., 2019), will be an important topic of ongoing study. The de- 
gree of imprinting may depend on the particular variants and the 
order in which they are introduced to the individual's immune 
system and the number of exposures, such as the number of 
vaccine doses received. Additional data for evaluating the 
magnitude of these effects and their consequences for protec- 
tion from infection are likely to become available in coming 
months, as individuals with different histories of SARS-CoV-2 
vaccination or viral variant infection become infected with the 
more highly mutated Omicron variant (https://covdb.stanford. 
edu/page/mutation-viewer/#omicron). As a practical consider- 
ation, the very-high spike-specific IgG concentrations generated 
by mRNA vaccination and periodic additional booster doses 
may be able to compensate for relatively decreased binding to 


(B) Representative in situ hybridization of an RNAScope control probe (left panels) and SARS-CoV-2 mRNA vaccine-specific probe (middle panels) within 
ipsilateral axillary core needle LN biopsies of female patients 7-60 days after second mRNA-1273 or BNT162b2 dose. Probe hybridization is indicated by red 
chromogen spots. IHC signal for spike antigen (right panels) is detected as granular brown staining. 

(C) Quantification of SARS-CoV-2 mRNA vaccine-specific probe-staining GCs in vaccinated LN biopsies. 

(D) Quantification of positive SARS-CoV-2 mRNA vaccine-specific probe spots per GC in vaccinee LNs. Error bars represent mean + SEM. 


(E) Spike-protein-positive GC quantification from IHC staining of vaccinee LNs. 

(F) IHC staining for spike (lower right panel) and nucleocapsid (upper panels and lower left panel) antigens in representative sections of COVID-19 patient 
peribronchial LNs. Nucleocapsid detection in primary (upper right panel) and secondary (upper left panel) LN follicles. 

(G) Due to the low frequency of detection of spike antigen in COVID-19 patient LNs, quantification is presented as the number of patients with positive staining in 
their LN specimens. 

(H) Quantification of the number of COVID-19 patients with LN follicles positive for nucleocapsid IHC staining. 

(l) Number and percentage of nucleocapsid-positive follicles by IHC in COVID-19 patient LNs. Error bars represent mean + SEM. 

(J) Spike concentration measured in plasma samples collected before and at several time points after BNT162b2 vaccination, with the red dotted line indicating 
the cutoff for positive. 

(K) Spike concentrations were measured in plasma samples collected from BNT162b2 vaccinees on DO (spike negative) or D28 (spike positive) spiked with 
different concentrations of recombinant spike protein. Black dotted line = cutoff for positive. 

(L) Spike concentration measured in plasma samples collected from BNT162b2 vaccinees on DO, D21, D22/23, and D28 mixed with the same plasma sample 
collected from one BNT162b2 vaccinee on D1 (spike positive). Black dotted line = cutoff for positive. 
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new viral variant antigens, potentially decreasing the public 
health impact of antibody response imprinting if vaccine boost- 
ing is widely adopted. 

We hypothesized that differences in the serological responses 
observed in SARS-CoV-2 infection compared with vaccination, 
particularly those related to variant-antigen-binding breadth, 
could be related to the anatomical sites where the viral antigens 
are encountered, the quantity of viral antigens in those anatom- 
ical sites, differences in the cell populations stimulated in 
secondary lymphoid tissues, and potential damage to immuno- 
logical tissues during infection. With CODEX multiplexed immu- 
nofluorescence microscopy and immunohistochemical micro- 
scopy, we identified follicular hyperplasia with robust axillary 
LN GOCs after mRNA (BNT162b2 or mRNA-1273) vaccination, 
containing CD21* follicular dendritic cell networks, BCL6* 
B cells, and PD-1* cells at significantly higher frequencies 
compared with those in peribronchial LNs of deceased COVID- 
19 patients. These findings demonstrate greater stimulation of 
GC B cells and Tfh cells in vaccination and normal functional or- 
ganization of GC follicular dendritic cells. Loss or impairment of 
GCs in patients with severe COVID-19 suggests that SARS- 
CoV-2 viral infection subverts the humoral immune response, 
by directly damaging immune cells or as a secondary effect of in- 
flammatory responses to infection (Feng et al., 2020; Kaneko 
et al., 2020). The observed extended presence of vaccine 
mRNA and spike protein in vaccinee LN GCs for up to 2 months 
after vaccination was in contrast to rare foci of viral spike protein 
in COVID-19 patient LNs. We hypothesize that the abundant 
spike antigen in the GCs of mRNA vaccine recipient LNs may 
contribute to the increased breadth of viral variant RBD binding 
by IgG seen after vaccination, potentially due to high antigen 
concentrations stimulating B cells with lower affinity for Wu- 
han-Hu-1 spike epitopes and better binding to variant epitopes. 
Persistent vaccine RNA and spike antigen at elevated concen- 
trations in vaccinee LNs could result in less strict selection for 
higher-affinity B cells in the immune response compared with 
situations where antigen is more limiting (Cirelli et al., 2019). 
However, our observation that all vaccine modalities (mRNA, 
adenoviral, and inactivated virus) stimulated greater viral variant 
breadth of IgG binding than infection could indicate that some 
other aspects of SARS-CoV-2 infection underlie these differ- 
ences, such as alteration of GC function. 

Prepandemic analysis of a model RNA vaccine for yellow fever 
virus in a rhesus macaque at 16 h postvaccination showed that 
vaccine RNA in LN cell suspensions was detected predomi- 
nantly in professional antigen-presenting cells including mono- 
cytes, classical dendritic cells, and B cells at this early time point 
(Lindsay et al., 2019). Data from follicular dendritic cells were not 
reported. Our histological data from SARS-CoV-2 mRNA-vacci- 
nated humans at considerably later time points (7-60 days post- 
second dose) show vaccine RNA almost entirely in GCs, distrib- 
uted primarily between the nuclei of GC cells, similar to the 
pattern seen by immunostaining for follicular dendritic cell pro- 
cesses or B cell cytoplasm. Additional co-localization studies 
with higher resolution may be required to determine more exactly 
which specific cell types harbor mRNA vaccine and spike anti- 
gen in humans following COVID-19 mRNA vaccination and infec- 
tion and may provide further mechanistic insights into the basis 
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for the differences in serological responses after vaccination 
compared with infection. 

At least some portion of spike antigen generated after admin- 
istration of BNT162b2 becomes distributed into the blood. We 
detected spike antigen in 9696 of vaccinees in plasma collected 
1-2 days after the prime injection, with antigen levels reaching as 
high as 174 pg/mL. The range of spike antigen concentrations in 
the blood of vaccinees at this early time point largely overlaps 
with the range of spike antigen concentrations reported in 
plasma in a study of acute infection (Ogata et al., 2020) although 
a small number of infected individuals had higher concentrations 
in the ng/mL range. At later time points after vaccination, the 
concentrations of spike antigen in blood quickly decrease 
although spike is still detectable in plasma in 6396 of vaccinees 
1 week after the first dose. A practical finding in our study is 
that the detection of spike antigen in plasma samples is impeded 
after second dose BNT162b2 vaccination, likely due to the for- 
mation of circulating immune complexes of anti-spike antibodies 
and spike protein, masking the antigen epitopes of the capture 
and detection antibodies that form the basis of antigen detection 
assays, similar to assay interference that has been reported for 
other diseases (Bollinger et al., 1992; Lima et al., 2014; Miles 
et al., 1993). 


Limitations of the study 
Data from SARS-CoV-2-infected clinical cohorts and vaccinated 
individuals in this study are observational. Longitudinal data for 
COVID-19 vaccine responses are derived predominantly from 
BNT162b2 mRNA vaccine recipients at Stanford, with data for 
the other four COVID-19 vaccines at a single postvaccination 
time point per individual. To make precise, internally controlled 
comparisons of polyclonal antibody responses to different viral 
variant antigens, we used multiplexed ECL assays of antibody 
binding to RBD, rather than virus neutralization assays, therefore 
our data do not reflect potentially functional antibodies binding 
to the spike N-terminal domain, or antibodies that may have 
other activities in vivo. In the analysis of imprinting of serological 
responses, plasma specimens were not available from the 
period after vaccination but before variant virus infection, pre- 
cluding direct comparison of antibody specificities pre- and 
post-infection. Additional epidemiological studies will be needed 
to evaluate the clinical impact of antibody response imprinting on 
susceptibility to infection by new viral variants and the severity of 
disease in infected patients. LN histological comparisons be- 
tween COVID-19 patients and vaccinees have the limitations 
that the infected patient specimens were limited to those with se- 
vere disease, the number of individuals analyzed was relatively 
low (seven COVID-19 patients and seven vaccinees), and the 
LN sampling was not done prospectively at predetermined 
time points after vaccination or infection. The serological anal- 
ysis in this study is of polyclonal antibody responses; evaluation 
of the clonal B cell and plasma-cell populations producing these 
antibody mixtures in comparable numbers of subjects in infec- 
tion and vaccination will likely be required for further mechanistic 
understanding. 

Taken together, these results underscore important differ- 
ences between SARS-CoV-2 antibody responses produced by 
vaccination versus infection. Key questions for the months and 
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years ahead include the duration of effective vaccine-stimulated 
serological responses after third-dose boosting or other 
repeated exposures, particularly for the recent Omicron variant 
and other variants that will emerge in the future, and the safety 
and efficacy of variant-targeting vaccine boosters in previously 
vaccinated or infected individuals. Further mechanistic investi- 
gations into the differences in antibody breadth elicited by vacci- 
nation and infection are needed to define the roles of T cell help, 
antibody affinity maturation, GC function, and innate immune re- 
sponses to vaccine components, as well as the cellular and sub- 
cellular distribution of vaccine RNA and expressed antigen in 
lymphoid tissues. Lessons from the antibody responses to the 
initial SARS-CoV-2 variants are likely to be important both for 
preparing for future additional variants of this virus, as well as 
improving vaccination strategies for other pathogens such as 
influenza virus. 
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Original code This paper; Zenodo https://doi.org/10.5281/zenodo.5854880 
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QuPath version 0.2.3 Bankhead et al. (2017) PMID: 29203879 
https://qupath.github.io/ 
Python version 3.7.10 Python Software Foundation https://www.python.org 
CODEX® MAV Akoya Biosciences https://help.codex.bio 
RESOURCE AVAILABILITY 


Lead contact 
Further information and requests for resources and reagents should be directed to the lead contact, Dr. Scott D. Boyd 
(publications_scott_boyd@stanford.edu). 


Materials availability 
This study did not generate new unique reagents. 


Data and code availability 

Raw data from all serology Figures have been deposited on Mendeley at https://doi.org/10.17632/j8r94pfrj6.1 and are publicly avail- 
able as of the date of publication. All original code has been deposited on the Zenodo platform at https://doi.org/10.5281/zenodo. 
5854880 and is publicly available at the Github repository https://github.com/boyd-lab/covid-infection-vs-vaccination as of the date 
of publication. Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon 
request. 


EXPERIMENTAL MODELS AND SUBJECT DETAILS 


Plasma and saliva samples from Stanford BNT162b2 vaccinees 

To study antibody responses after first, second, and third dose vaccination with BNT162b2, we collected longitudinal blood and 
saliva samples from 59 vaccinees (29 were women, 27 were men, 3 were unknown, and all donors were adults between the ages 
of 19 to 79 years). Baseline blood samples were collected on day 0 before or immediately after the first vaccine dose. Individuals 
received their second dose on day 21 and a third dose about 9 months after the prime. Blood sample collection after prime vacci- 
nation was scheduled for days 1, 7, 21, 22, 28, and 42 with blood draws + one day from the assigned time point, or days 90, 120, or 
210 + 1 week from the assigned time point. In addition, blood samples were collected on days 0 to 3, 7 to 10, 21, and 28 after the third 
vaccine dose. Saliva samples were collected on day 42 after the first dose, as well as before, and 1 or 2, 3 or 4, 7, and 21 days after the 
third vaccine dose. Peripheral blood was collected in vacutainer cell preparation tubes (CPTs) containing sodium citrate. Plasma was 
isolated and stored at -80°C. Saliva was collected from study participants, centrifuged, and supernatants were stored at —80°C. All 
BNT162b2 vaccine study participants provided informed consent under Stanford University Institutional Review Board approved 
protocol IRB-55689. 


Plasma and saliva samples from Stanford COVID-19 patients 

Blood and saliva samples were collected between March and December 2020 from COVID-19 patients who reported to Stanford 
Healthcare-associated clinical sites. SARS-CoV-2 infection was confirmed for all patients by RT-qPCR of nasopharyngeal swabs 
as described (Corman et al., 2020; Hogan et al., 2020). Blood samples were collected in heparin- or EDTA-coated vacutainers. After 
centrifugation for collection of plasma or saliva, samples were stored at —80°C. The use of these samples for antibody testing was 
approved by the Stanford University Institutional Review Board (Protocols IRB-48973 and IRB-55689). 

Stanford COVID-19 patient cohort 1 included 530 plasma samples collected from 100 moderately to severely ill COVID-19 patients 
(52 were women and 48 were men; ages ranged from 1 to 95 years; 23 were outpatients, 33 were admitted to hospital without needing 
ICU care, 20 were treated in the ICU and 24 died of COVID-19) between March 2020 and August 2020. ELISA serology data for Wu- 
han-Hu-1 SARS-CoV-2-specific antibodies in these specimens have been reported previously (Róltgen et al., 2020). 

Stanford COVID-19 patient cohort 2 was included as a validation cohort of 87 samples from 74 mostly mildly ill patients who had 
blood sample draws between March and December 2020, at approximately 21 days (n = 15 samples), 1 month (n = 23 samples), 
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3 months (n = 27 samples) and 7 months (n = 22 samples) after positive RT-qPCR testing for SARS-CoV-2 infection. Of those patients, 
37 were women, 34 were men and 3 were unknown. Donors were 19 to 72 years of age and in terms of disease severity, 59 were 
mildly ill, 6 were moderate ill and 9 had a severe/critical disease course. Specimen time points were selected to match those of Stan- 
ford BNT162b2 vaccinee sample collections. Saliva samples were collected from five COVID-19 patients. 

Stanford SARS-CoV-2 variant infection cohort blood samples were collected from COVID-19 patients during acute infection with 
SARS-CoV-2 Alpha (n = 7) or Delta (n = 34) variants. Samples were from 20 women and 21 men, all between 2 and 92 years of age. 
SARS-CoV-2 genotyping data were obtained using a multiplex, mutation-specific RT-qPCR targeting N501Y, E484K, and L452R, as 
previously described (Wang et al., 2021). Samples from the first multiplexed reaction suspected to contain the Alpha variant were 
analyzed with a second confirmatory genotyping RT-qPCR assay to detect mutations encoding the N501Y amino acid change, 
as described (Dashdorj et al., 20218). 


Plasma samples from Mongolian vaccinees 

To study SARS-CoV-2 variant-specific IgG responses elicited by different COVID-19 vaccines, we tested plasma samples collected 
in July 2021 from 196 Mongolian vaccine study participants (109 were women, 87 were men, all were adults between 20 and 85 years 
of age) who had been fully vaccinated with one of four COVID-19 vaccines: BioNTech-Pfizer BNT162b2 (n = 47), AstraZeneca ChA- 
dOx1-S (n = 50), Sputnik V Gam-COVID-Vac (п = 45) and Sinopharm BBIBP-CorV (n = 54). Participants were recruited by public 
announcement and volunteers were enrolled after signing the consent form approved by the Ethics Review Board at the Ministry 
of Health of Mongolia. SARS-CoV-2 pseudotyped virus neutralization and RBD-ACE2 blocking data on the same samples have 
been reported previously (Dashdorj et al., 20218). Peripheral blood was collected іп CPT, centrifuged for collection of plasma, 
and stored at —80°C. 


Healthy human control (HHC) plasma and saliva 

37 plasma and 20 saliva samples from HHCs collected before the onset of the COVID-19 pandemic for studies at the Sean N. Parker 
Center for Allergy & Asthma Research were used to verify pre-pandemic antibody binding concentrations to the different coronavirus 
antigens, and manufacturer-provided cutoffs for positive serology assay results. Use of these samples was approved by the Stanford 
University Institutional Review Board (Protocols IRB-8629 and IRB-60171). No demographic information was available for these 
samples. 


Axillary LN core biopsies and post-mortem peribronchial LN tissues 

To analyze and compare GC architecture in response to COVID-19 vaccination and SARS-CoV-2 infection, we collected axillary LN 
core needle biopsies from BNT162b2 or mRNA-1273 vaccinees, and excised post-mortem peribronchial LNs from patients who died 
of COVID-19. For the selection of vaccinee tissues, we performed a retrospective search of our pathology archives and medical re- 
cords between January 2021 and June 2021 for female patients who received either mRNA-1273 or BNT162b2 vaccination and sub- 
sequently underwent an ipsilateral axillary LN core needle biopsy due to mammographic findings and routine clinical care. Seven 
patients underwent biopsy one to eight weeks after vaccination with their second dose of mRNA vaccine. Three unvaccinated fe- 
males undergoing axillary LN core biopsy for routine clinical care and mammographic findings served as controls. We included seven 
peribronchial LNs from two female and five male patients who died of COVID-19 before August 2020, one to three weeks after symp- 
tom onset. Control post-mortem peribronchial LN biopsies were from pre-pandemic patients who died of non-COVID-19 causes. 
Autopsies were done by the Arizona Study of Aging and Neurodegenerative Disorders Brain and Body Donation Program (Beach 
et al., 2015). Analysis of these tissues was approved by Stanford University Institutional Review Board Protocol IRB-48973. 


METHOD DETAILS 


MSD ECL binding assays 

Plasma samples from vaccinees and COVID-19 patients were heat-inactivated at 56°C for 30 minutes and tested using multiplexed 
ECL detection in a 96-well plate format with MSD® V-PLEX® serology panels and instrumentation according to the manufacturer's 
instructions. V-PLEX COVID-19 Coronavirus Panel 2 kits were used to detect IgM, IgG, and IgA antibodies to SARS-CoV-2 N, S1 
NTD, RBD, and spike antigens and to spike proteins of SARS-CoV-1 and other HCoVs including HCoV-OC43, HCoV-HKU1, 
HCoV-NL63, and HCoV-229E. V-PLEX SARS-CoV-2 Panel 9 and 11 kits were used to determine IgG antibody concentrations 
and RBD-ACE2 blocking antibody percentages to different SARS-CoV-2 variant RBDs, with Alpha, Beta, Gamma, Epsilon, Kappa, 
Eta/lota, B.1.526.2, P.3 and Wuhan-Hu-1 present in both panel 9 and 11, and with B.1.214.2 in panel 9 and Delta in panel 11. V-PLEX 
SARS-CoV-2 Panel 20 kits were used to determine IgG antibody concentrations to Alpha, Beta, Gamma, Delta, and Wuhan-Hu-1 
SARS-CoV-2 variant spike proteins. Plasma samples were analyzed in duplicate at a 1:5,000 (for IgG binding assays) or a 1:10 
(for RBD-ACE2 blocking assays) dilution in MSD diluent. Coronavirus-specific antibodies were detected with anti-human IgM, 
IgG, or IgA antibodies, or indirectly with human ACE2 protein (for RBD-ACE2 blocking assays) conjugated to SULFO-TAG™ ECL 
labels and read with a MESO® QuickPlex® SQ 120 instrument. Cutoff values for positive antibody test results for Wuhan-Hu-1 
antigens were determined by the manufacturer based on sera from 200 pre-pandemic healthy adults and 214 PCR-confirmed 
COVID-19 patients. We tested an additional 37 healthy adult pre-pandemic plasma specimens to evaluate the manufacturer's cutoff 
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values, and to determine cutoffs for positive binding to variant virus antigens, defined as the mean plus three standard deviations of 
the results from the pre-pandemic specimens. Antibody binding ratios for Wuhan-Hu-1 and viral variant antigens were only calcu- 
lated for specimens that were above the cutoff values for positive results. Saliva samples were analyzed in duplicate at a 1:5 dilution 
in MSD diluent 2. Each plate contained duplicates of a 7-point calibration curve with serial dilution of a reference standard, a blank 
well and three positive control samples. Calibration curves were used to calculate antibody unit concentrations by backfitting ECL 
signals measured for each sample to the curve. 


MSD ECL spike antigen detection 

SARS-CoV-2 spike antigen was quantified in plasma samples using an antigen capture ECL immunoassay platform (Meso Scale Dis- 
covery). S-PLEX& SARS-CoV-2 spike Kit assays were performed according to manufacturer instructions. A 7-point calibration curve 
and negative control consisting of assay diluent were run in duplicate on each plate. Plates were read using a MESO QuickPlex SQ 
120 instrument. Raw signals were converted to a concentration based on linear regression to the 7-point calibration curve. Recom- 
binant SARS-CoV-2 spike protein used for plasma spiking experiments was made by ATUM (https://www.atum.bio/). 


Histology, immunohistochemistry, and in situ hybridization 
LN core needle and autopsy tissue samples were fixed in formalin and embedded in paraffin (FFPE), and sectioned. Once unstained 
slides were generated and initial Hematoxylin and Eosin (H&E) stained sections were analyzed, two distinct 0.6 mm areas from each 
LN sample were cored out of each tissue block and re-embedded to construct a tissue microarray (TMA). Immunohistochemistry was 
performed on four-micron sections using standard automated or manual methods including deparaffinization, peroxidase blocking, 
antigen retrieval, primary and secondary antibody incubation, detection with 3,3'-Diaminobenzidine (DAB) development, and coun- 
terstaining. Assays were performed on Roche Ventana (Tucson, AZ) Ultra instruments using Ventana Optiview detection, or Leica 
(Buffalo Grove, IL) Bond III instruments using Leica Polymer Refine detection or manually using Dako (Carpenteria, CA) Target 
Retrieval (TR) and Liquid DAB+ Substrate Chromogen System with ImmPress (Vector, Burlingame, CA) secondary antibodies. 
For in situ hybridization, manual methods were used as previously described (Cloutier et al., 2021), using manufacturer-recommen- 
ded protocols with the RNAScope 2.5 HD Assay-RED kit and probes from Advanced Cell Diagnostics (Newark, CA). Two SARS- 
CoV-2 vaccine probes were developed to target bases 101-1143 of the spike encoding sequence of the BNT162b2 vaccine or bases 
101-1488 of the spike encoding sequence of the mRNA-1273 vaccine. Both probes recognized SARS-CoV-2 mRNA vaccine, thus 
only the SARS-CoV-2 vaccine probe recognizing bases 101-1488 of mRNA-1273 vaccine are presented. To assess the specificity of 
SARS-CoV-2 RNAScope vaccine probes, they were tested against SARS-CoV-2 infected placental tissue, in addition to staining for 
SARS-CoV-2 viral probe which targets bases 21631-23303 of the S-gene. SARS-CoV-2 RNAScope vaccine probes did not recog- 
nize SARS-CoV-2 virus. In addition, for each tissue tested there were internal negative control areas which did not react with SARS- 
CoV-2 vaccine probe. 


Whole slide imaging (WSI) and quantitative image analysis 

Whole slide images of immunohistochemistry-stained slides for BCL6, PD-1 and CD21 were scanned at 40X magnification (0.25 uM 
per pixel) on an Aperio AT2 scanner (Leica Biosystems, Nussloch, Germany) in ScanScope Virtual Slide (SVS) format. Total lymphoid 
tissue and individual GCs were annotated as regions of interests (ROIs) using the QuPath open-source WSI software (Bankhead et al., 
2017). GCs were defined as B cell areas with CD21* follicular dendritic cell networks and BCL6+ nuclei. Primary follicles were defined 
as B cell areas with CD21* follicular dendritic cell networks without BCL6-- nuclei. For each GC ROI, the standard positive cell detec- 
tion function was used to identify positive and negative cells with a single threshold of 0.2 when scoring the cell DAB OD mean. For 
autopsy BCL6 analysis (COVID-19 and control) the threshold was adjusted to 0.05 given the dim expression of BCL6. For each ROI, 
the area in mm?, number positive cells per mm?, and percent positive cells were calculated by QuPath. 

ATMA slide containing two distinct 0.6 mm cores of each of the mRNA vaccine (n = 7) and vaccine control (n = 3) biopsies hybrid- 
ized with SARS-CoV-2 mRNA RNAScope vaccine probe was scanned at 40x magnification (0.25 uM per pixel) on an Aperio AT2 
scanner in SVS format. Total lymphoid tissue and individual GCs were annotated as ROIs using QuPath. For each GC ROI, the num- 
ber of spots/clusters of RNAScope probe was detected using the QuPath subcellular detection option per manufacture instructions. 
Our detection parameters were as follows: detection threshold = 0.6; split by intensity. Our split and cluster parameters were as fol- 
lows: expected spot size = 2 um?, min spot size = 2 um?, and max spot size = 3 um?. 


Co-detection by indexing (Codex) 

All antibodies used for Codex were first screened, titrated, and validated by individual staining on FFPE human tonsil tissue samples. 
Standard manual immunohistochemistry was used to crossvalidate antibodies with the same, nonconjugated antibody clones listed 
in the key resource table. Tissue preparation for Codex was undertaken by obtaining 8-ium-thick sections from FFPE tissue blocks 
that were immobilized on charged square glass coverslips coated with polylysine (Electron Microscopy Sciences Hatfield, PA) pre- 
pared according to the manufacturer's instructions. The coated glass coverslips were stained with a cocktail containing nucleotide- 
barcoded primary antibodies. The coverslips underwent nuclear staining (DAPI) and were loaded on the stage of an automated 
inverted fluorescence microscope connected to the robotic fluidic system known as CODEX (Akoya Biosciences, Marlborough, 
MA). In addition to the nuclear stain, fluorophore-tagged complementary nucleotide sequences (reporters) were used to iteratively 
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reveal three antibodies at a time per cycle. Two additional blank cycles at the beginning and at the end of the antibody reveal cycles 
were added for purposes of subtraction of autofluorescence background. Automated image acquisition, processing, segmentation, 
and fluidics exchange were conducted using an Akoya Codex instrument and Codex driver software (Akoya Biosciences) (Black 
et al., 2021; Goltsev et al., 2018; Schürch et al., 2020). Composite images providing multiplex fluorescent signals of stained cell types, 
cellular niches, and tissue architecture were captured using CODEX® MAV and FIJI software (Akoya Biosciences). 


QUANTIFICATION AND STATISTICAL ANALYSIS 


Statistical tests were performed in R using base packages for statistical analysis and the ggplot2 package for graphics. Box-whisker 
plots show median (horizontal line), interquartile range (box), and the end of the lower whisker representing the smallest observation 
greater than or equal to the 2596 quantile minus 1.5 times the interquartile range, and the end of the upper whisker representing the 
largest observation less than or equal to the 7596 quantile plus 1.5 times the interquartile range. In serological analyses where sta- 
tistical significance was tested, significance was defined as: ***p value < 0.001; **p value < 0.01; *p value < 0.05. 

Forthe principal component analysis (PCA) we log-transformed, calculated z-scores, and ran PCA on MSD antibody concentration 
measurements or Wuhan-Hu-1/variant RBD IgG concentration ratios from a reference time point after COVID-19 vaccination or 
SARS-CoV-2 infection using Python v3.7.10 and packages numpy v1.19.1, pandas v1.2.5, and scikit-learn v1.0. We then applied 
these transformations to matching data from all other time points, enabling us to visualize the change over time in these serology 
measurements on a consistent PCA reference. Plots were created with Python packages matplotlib version 3.3.2 and seaborn 
version 0.11.2. 

For the analysis of the homogeneity or dispersion of serology measurements in groups differing by vaccination or infection status 
(Figures 2C and 2D) at a particular time point, we plotted each group's distribution of Euclidean distances to its centroid (calculated 
with Python package scipy version 1.6.2). These distance distributions were consistent when calculated in the raw measurement 
data space of arbitrary units (AU) for the MSD ECL assay or in the transformed PC1 and PC2 space after embedding the raw mea- 
surements into the PCA space created from one reference time point. 

Ratios of concentrations of IgG binding to Wuhan-Hu-1 RBD compared to variant virus RBDs were plotted for specimens with IgG 
binding above the cutoff for positive binding to Wuhan-Hu-1 RBD (Figures 3, 4 and S4) to avoid distortion of ratios by samples without 
specific binding. Ratios of IgG binding to spike antigens were calculated in a similar manner (Figure S4C, lower panel). Corresponding 
IgG concentrations for samples used to calculate ratios were plotted for reference (Figures 3A, 4A, S4A and S4C upper panel). To 
quantify serological imprinting from prior Wuhan-Hu-1 antigen exposure from vaccination on subsequent responses to breakthrough 
infection with the Delta variant, we first computed the ratio of Wuhan-Hu-1 RBD-binding level to Delta RBD-binding level. Here, a ratio 
of one indicates even preference, while ratios greater than one indicate preferential binding of Wuhan-Hu-1 over Delta. Each ratio is 
symmetric with its inverse; for example, a Wuhan-Hu-1/Delta binding ratio of 4/5 indicates the same degree of preference for Delta 
binding as the ratio 5/4 indicates for Wuhan-Hu-1 binding preference. We then log-transformed the ratios, which sets the even pref- 
erence level at zero, with positive values corresponding to Wuhan-Hu-1 preference, and makes these values symmetric around zero 
(e.g., a value of -0.2 indicates the same level of preference for Delta as +0.2 does for Wuhan-Hu-1 binding). Finally, we rescaled the 
negative and positive values separately to the ranges -100% to 0 and 0 to +100%. The resulting magnitudes are binding preferences 
relative to the maximum binding preference observed for a particular variant, including data from individuals only exposed to Wuhan- 
Hu-1 or individual variants. In particular, -10096 refers to the maximum observed binding preference towards Delta. We plotted the 
distributions of these binding preference levels for BNT162b2 vaccinees without known prior SARS-CoV-2 infection, for Delta infec- 
tion cases with no recorded prior infection or vaccine exposure, and for Delta breakthrough infections following Wuhan-Hu-1 wild- 
type vaccination. An example binding preference level of +30% for a Delta breakthrough infection case suggests this individual is at 
3096 of the most imprinted state for Wuhan-Hu-1 preference. A binding preference of -2096 would instead suggest that, following 
infection, this individual lost imprinting to Wuhan-Hu-1 and gained preference for binding the Delta variant. 
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Figure S1. Anti-SARS-CoV-2 Ig antibody responses in plasma and saliva following BNT162b2 vaccination, related to Figure 1 
Anti-SARS-CoV-2 N, RBD, and spike (S) IgG (A), IgM (B), and IgA (C) responses are shown for plasma from individuals who received BNT162b2 prime (DO, n = 59) 
and second dose (D21, n = 58) vaccination. Box-whisker plots of the WHO binding arbitrary unit (BAU/mL) anti-SARS-CoV-2 concentrations show the inter- 
quartile range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 
quantile. Comparisons between groups of previously SARS-CoV-2-infected (CoV-2+) versus noninfected individuals, and female versus male were by the two- 
sided Wilcoxon rank sum test; comparison between age groups («40; 40-60; >60 years) was done using pairwise Wilcoxon rank sum test with Bonferroni 
correction. *p < 0.05, **р < 0.01, ***р < 0.001. 

(D) Anti-SARS-CoV-2 N, RBD, and S IgG concentrations in BAU/mL are shown for saliva from individuals who received BNT162b2 prime/boost and third dose 
vaccination (upper left panel). Anti-SARS-CoV-2 N, RBD, and S (upper right panel) concentrations in BAU/mL, as well as anti-SARS-CoV-1 and anti-HCoV-OC43, 
-HKU1, -NL63, and -229E S IgG (lower panel) concentrations in MSD arbitrary units (AU/mL), are shown for saliva collected on D42 after BNT162b2 prime 
vaccination (vaccinee), around D42 post-symptom onset for COVID-19 patients (CoV-2+), and before the onset of the COVID-19 pandemic for pre-pandemic 
healthy human controls (Pre-pan). Box-whisker plots of anti-SARS-CoV-2 IgG concentrations show the interquartile range as the box and the whisker ends as the 
most extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 quantile. Statistical test for significance between groups 
(CoV-2+; Pre-pan, vaccinee) was performed using pairwise Wilcoxon rank sum test with Bonferroni correction. *p « 0.05, **p « 0.01, ***р < 0.001. 
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Figure S2. The magnitude of antibody responses is not correlated with reported vaccine-associated side effects (SEs), related to Figure 1 
(A) Frequency of site-specific and systemic vaccine-associated SEs after prime (light green) and second dose (dark green) BNT162b2 vaccination. 

(B) Box-whisker plots of the MSD AU/mL anti-SARS-CoV-2 IgG concentrations in BNT162b2 vaccinee plasma collected on D28 postvaccination show the in- 
terquartile range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 25% quantile and above the 75% 
quantile. For a given SE (rows), vaccinees were grouped according to no SE reported (“No,” colored in blue) or SE reported ("Yes," colored in orange). Vaccinees 
where SEs were unknown are shown as white boxplots. 
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Figure S3. BNT162b2 vaccination produces less broad serological responses to endemic human coronaviruses (HCoVs) compared with 
SARS-CoV-2 infection, related to Figure 2 

(A and B) Anti-SARS-CoV-1 spike and anti-HCoV-OC43, -HKU1, -NL63, and -229E spike IgM, IgG, and IgA antibody responses are shown for individuals who 
received BNT162b2 prime (DO) and boost (D21) vaccination doses and for COVID-19 patients. (A) The heatmap shows the development of antibody responses in 
longitudinal samples from vaccinees/patients collected at/during DO, D7/week 1, D21/weeks 2 and 3, D28/week 4, D42/weeks 5 and 6, and D90/> week 7 after 
vaccination/COVID-19 symptom onset (x axis). The color scale encodes the median values of одо MSD AU/mL concentrations. (B) Box-whisker plots show the 
development of antibody responses in longitudinal samples from vaccinees/patients collected at/during DO, D7/week 1, D21/weeks 2 and 3, D28/week 4, D42/ 
weeks 5 and 6, and D90/> week 7 after vaccination/COVID-19 symptom onset (x axis). Box-whisker plots show the interquartile range as the box and the whisker 
ends as the most extreme values within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. Statistical test: pairwise Wilcoxon 
rank sum test with Bonferroni correction. *p « 0.05, **p « 0.01, ***p « 0.001. Individuals were classified as vaccinees who have not been previously exposed to 
SARS-CoV-2 (vaccinees), outpatients (Outpt) and hospital-admitted patients (Admit), and ICU patients and those who died from their illness (Death). 

(C) Box-whisker plots show anti-SARS-CoV-1 spike and anti-HCoV spike antibody responses in plasma samples from individuals who received BNT162b2 prime 
(DO, n = 59 individuals), second dose (D21, n = 58 individuals), and third dose (around month 9, n = 36 individuals) vaccination. Box-whisker plots show the 
interquartile range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 25% quantile and above the 75% 
quantile. 
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Figure S4. Greater breadth of IgG binding to SARS-CoV-2 variant RBDs following BNT162b2 vaccination compared with infection with 
Wuhan-Hu-1 SARS-CoV-2 (validation cohort), related to Figure 3 

(A and B) Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG responses are shown for Stanford individuals who received BNT162b2 vaccination and for 
Wuhan-Hu-1-infected COVID-19 Stanford patient cohort 2 at different time points after vaccination/COVID-19 symptom onset. Box-whisker plots show the 
interquartile range as the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 2596 quantile and above the 7596 
quantile. Significance between groups was tested with two-sided Wilcoxon rank sum test. *p < 0.05, **р < 0.01, ***p < 0.001. 

(A) Anti-RBD IgG concentrations. 

(B) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration. 

(C) Anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant spike IgG responses as anti-spike IgG concentrations (upper panels) and as ratios of anti-Wuhan-Hu-1 to 
variant spike IgG concentration (lower panels) are shown for Stanford individuals who received BNT162b2 vaccination and for Wuhan-Hu-1-infected COVID-19 
Stanford patient cohorts 1 and 2 samples. Box-whisker plots show the interquartile range as the box and the whisker ends as the most extreme values within 1.5 
times the interquartile range below the 2596 quantile and above the 7596 quantile. Significance between groups was tested with two-sided Wilcoxon rank sum 
test. *p < 0.05, **р < 0.01, ***р < 0.001. 

(D) Percentage blocking of ACE2 binding to RBD of specified viral variants by plasma antibodies of BNT162b2 vaccinees and Stanford patient cohort 2 samples. 
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Figure S5. Anti-SARS-CoV-2 RBD IgG signatures following BNT162b2 vaccination and SARS-CoV-2 infection, related to Figure 5 

(A) Ratios of anti-Wuhan-Hu-1 to variant RBD IgG concentration are shown for Stanford individuals who received BNT162b2 vaccination at different time points 
after the second dose (D21, n = 58 individuals) and third dose (around month 9, n = 36 individuals) vaccination. Box-whisker plots show the interquartile range as 
the box and the whisker ends as the most extreme values within 1.5 times the interquartile range below the 25% quantile and above the 75% quantile. 

(B) Principal component analysis (PCA) of anti-SARS-CoV-2 Wuhan-Hu-1 and viral variant RBD IgG concentrations across Stanford BNT162b2 vaccinees, 
Stanford COVID-19 patient cohort 2, and SARS-CoV-2-variant-infected patients. 
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Figure S6. Disrupted LN GCs in COVID-19 patients versus mRNA vaccinees, related to Figure 6 

(A) LN GC histology for COVID-19 patients (left) and mRNA vaccinees (right) evaluated with four-color codetection by indexing (Codex) immunofluorescence 
analysis for CD20 (red), CD3 (blue), BCL6 (magenta), and CD21 (yellow) markers of B cells, T cells, GC B cells (or T follicular helper cells), and follicular dendritic 
cells, respectively. 

(B) Representative CD21 immunohistochemistry of secondary (left) and primary (right) follicles of four autopsy patients who died of COVID-19 and two control 
autopsy patients. 


